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Abstract 

Geometric phases of simple harmonic osciUator system are studied. Complete 
sets of " eigenfunctions" are constructed, which depend on the way of choosing clas- 
sical solutions. For an eigenfunction, two different motions of the probability dis- 
tribution function (pulsation of the width and oscillation of the center) contribute 
to the geometric phase which can be given in terms of the parameters of classical 
solutions. The geometric phase for a general wave function is also given. If a wave 
function has a parity under the inversion of space coordinate, then the geometric 
phase can be defined under the evolution of half of the period of classical motions. 
For the driven case, geometric phases are given in terms of Fourier coefficients of 
the external force. The oscillator systems whose classical equation of motion is 
Mathieu's equation are perturbatively studied, and the first term of nonvanishing 
geometric phase is calculated. 
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1. Introduction 

The harmonic oscillator is a system whose wave functions can be given in terms of so- 
lutions of classical equation of motion. This fact has been recognized by Lewis fll who 
found, in an application of asymptotic theory of Kruskal , that there exists a quantum 
mechanically invariant operator which has been used to construct wave functions of (gen- 
eralized) harmonic oscillator systems. An alternative and simple (at least, conceptually) 
way to construct a complete set of "eigenfunctions" in terms of classical solutions is to use 
the Feynman path integral method 0. Based on these progresses, the geometric phases 
0, Q for the eigenfunctions of time-periodic harmonic oscillator system has been studied 
ip. The simple harmonic oscillator (SHO) of constant mass and constant frequency is 
a system where any real number can be a period of the Hamiltonian and two linearly 
independent classical solutions are finite all over the time. In view of the general analyses 
P], for SHO there is a representation where the geometric phase can be defined for an 
evolution of any real number period. For the driven case the phases can not be defined if 
a particular solution diverges 0, . 

We note that the geometric phase of SHO has been studied in In considering a 

spin-1/2 particle in a harmonic potential and a nonuniform magnetic field, Pati and Josh 



consider a wave function of SHO whose center oscillates [10] and calculate the geometric 



phase. By adopting the fact that a Gaussian wave packet could be a wave function of a 



harmonic oscillator system, geometric phase has been calculated [O for the wave packet 



whose width oscillates and relation between the phase and Hannay's angle []T2|] is given 
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for the wave packet. 

In this paper, the geometric phases in SHO system will be studied systematically. 
Making use of the general formalism 0, we will construct a complete set of eigenfunc- 
tions and show that, for an eigenfunction, the geometric phase comes from two different 
motions of the probability distribution function: the oscillation of the center and the 
pulsation of the width. The driven case will be also considered. For a harmonic oscillator 
system whose classical equation of motion is Mathieu's equation, the geometric phase 
for an eigenstate will be perturbatively studied in a representation which reduces to the 
"stationary representation" of SHO in a limit. 

In the next section, by applying the general formulas 0, the geometric phases for 
the eigenstates of SHO with driving force will given by assuming the existence of periodic 
classical solution. In section 3, the geometric phase of an eigenfunction will be interpreted 
as the contributions from the two different motions. In section 4, geometric phase will 
be calculated for a general wave function of SHO. In section 5, it will be shown that, if a 
wave function of SHO has a parity under the inversion of space coordinate, the geometric 
phase of the wave function can be defined under the evolution of half of the period of 
classical motions in SHO. In section 6, the geometric phase for an eigenstate in Mathieu's 
oscillator system will be perturbatively treated. 
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2. General formula 

Classical equation of motion of the SHO with driving force F{t) is written as 

x + w'x = ^. (1) 

If the two linearly independent homogeneous solutions and a particular solution of (1) 
are denoted as u(t),v{t) and Xp(t), respectively, by introducing 6(t) defined through the 
relation 

S = ^Mw^Xp - ^Mxp, (2) 



the wave functions satisfying Schrodinger equation can be written as P, [13 

i^mix, t) = < x\m > 

where p{t) and are defined as 



p{t) = ^u\t)+v^{t) Q = M{uv - vii). (4) 

With given u{t),v{t) and Xp{t), the wave functions of (3) satisfying orthonomality condi- 
tions 

/oo 
^l^mdx = Sn,m, (5) 
-oo 

form a complete set as can be seen from ^ 

J2'^m{Xb,'tb)i'^i^a,'ta) ^ S{Xb- Xa) aS U i ta- (6) 



Considering this completeness, ipmixjt) will be referred to as an eigenfunction. 
The geometric phase for the eigenfunction iljm{x,t) is written as p 

l^ir') = -(m + -)[ (^-—^ + ^^p^)dt+—i xldt 



where r' is a period needed for the definition of geometric phase. 

Without losing generality, the two linearly independent homogeneous solutions u{t),v{t) 
can be written as 



u{t) = cos wt v{t) = Csm{wt + p), (8) 

where C is a nonzero real number and a real number (3 is not one of {2n + l)7r/2 (n = 
0, ±1, ±2, ■ ■ ■). Then Q is simply written as 

n = MwC cos (3. (9) 

For the study of geometric phase, we only consider periodic F satisfying F{t + Tf) = F{t). 
For the existence of geometric phase, the Xp should be finite all over the time which will 
be assumed from now on. If there exist two positive integers N,p of no common divisor 
except 1 such that tq/tj = p/N, by defining tq = 2ti/w, Ntq can be the r' in (7). The 
periodic F{t) can be written as 



oo 



Fit)= Yl (10) 
n=— oo 

where 

U = - r F{t)e''^^f'dt (11) 



Tf JO 



with Wf = 27r/r/. The Xp{t) finite all over the time can be written as 

oo f 

where D and (f are real numbers, li p — 1, must be zero for the finiteness of Xp. The 
minimum period of Xp is Ntq (= pTf) in general so that 

Xp{t + NTo)^Xp{t). (13) 

After some algebra, from (7) one may find the geometric phase of the eigenfunction 
ipn{x,t) under A'"ro-evolution as 

F/.r N / l^,,,l-2Ccosf3 + C\ 
= nin + -)N[ ^-^ ] 

It should be mentioned that different choice of (C, /3, D, cp) gives different representation 
(set of eigenfunctions) , and thus different geometric phases. In the next section, an 
explanation for the fact that the geometric phase does not depend on (p will be given. 

There exist the representations where geometric phase can not be defined even with 
periodic Xp. For example, in the representations of D 7^ 1, if t//to is an irrational number 
geometric phases can not be defined under any evolution of finite time. However, for the 
SHO with driving force, there is a special representation of C = 1, /? = and D = where 
geometric phase can be defined for any periodic Xp. In this representation, geometric phase 
under r/-evolution is written as 



For the SHO without driving force, if we take C = 1 (3 = and Xp = 0, the probabihty 
distributions of the eigenstates are static, and thus the representation of such choice of 



parameters will be referred to as the stationary representation |T3|. 

In the representation of C = 1, (3 = the geometric phase for the driven harmonic 
oscillator is considered in |P, but the given results there do not agree with (14) and (15). 

3. Fictitious particular solution and geometric phase 

From now on we only consider case of F = 0. 

If we take Xp = 0, the geometric phase for ipn under the ro-evolution is given as 

, 1,^ 1 -2C cos/? + 
7n(ro) = vr(n + -)[ ]. (16) 

Since § 

roo 

'?/'*(x, t)xipn{x, t)dx = Xp{t), (17) 



if we take Xp = 0, the center of the probability distribution function of an eigenfunction 
does not move. In the stationary representation, the probability distribution given from 
an eigenstate does not move as time passes, so that the geometric phase for an eigenstate 
which can be defined for evolution of any period is zero |14]. It has been noted |l^ 



that, except for the stationary representation, the widths of probability distributions of 
the eigenstates pulsate as time passes by, which may be the physical origin of the geometric 
phases given in (16). The geometric phase for a Gaussian wave packet whose center of 



probability distribution does not move has been considered in [11|. In order to compare 
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with the the result of Ge and Child |]TT[, by defining 

2Ti p 



i ro a , l-2Ccosf3 + C^ 
-dt = vr 



one may find the relation 



2 Jo a + a* 2C cos j3 

which gives the geometric phase 7o(to) in agreement with the result in (16). 

Even for the SHO without driving force, one can take the fictitious particular solution 

as 

Xp = D co?,{wt + if). (19) 

It may be straightforward, from (3), to find the eigenfunctions for SHO with the Xp in 
(19) and check that the eigenfunctions indeed satisfy the Schrodinger equation. To make 
the point clear, we only explicitly written down the eigenfunctions of C = 1 and /? = 



ijjn = J _ ^ exp[ial[-xD sm{wt + ^p) + \d'^ sin(2wt + y?)] - i{n + \)wt] 

X exp[ ^]if„(ao(a; - D cos{wt + if ))), (20) 



where ao = y Mw/h. From the general formula of (14), one can find that, the geometric 
phase for an eigenfunction ipn in (20) under the ro-evolution is given as nMwD"^ /h (= 
vrag-D^). ipo in (20) coincides with the wave function whose geometric phase is evaluated 



in the study of a spin-1/2 particle in a nonuniform magnetic field [^. The fact that the 
geometric phase does not depend on if is simply understood from that (f determines the 
initial point on a cycle while the geometric phase is calculated through the integral over 
a cycle. 



This section may be summarized as follows; For an eigenstate in (3) of SHO, with the 
choice of solutions of classical equation of motion in (8,19), the geometric phase under 
To-evolution is written as 

T„(T.)=.(»+^)i — — i+-ir^- Pi) 

The first term in the right-hand side of (21) is from the pulsation of width of the probabihty 
distribution of the eigenstate and the second term is from the oscillation of the center of 
the probability distribution. 

4. Geometric phase for a general wave function of 
SHO 

Due to the completeness, any wave function satisfying Schrodinger equation of (driven) 
SHO may be written as a linear combination of eigenfunctions with a specific choice of 
classical solutions: 

oo 

^^J2BnM^,t). (22) 
n=0 

For simplicity, we only consider SHO {F — 0) system. The classical solutions will be 
chosen as in (8,19). Then one may easily find that every eigenfunction satisfies the 
relation ipniXjTo + t) = exp[— i(2n + l)7i]il!n{x,t). Therefore, under the ro-evolution, the 
geometric phase can be defined for any wave function of SHO. 

For the explicit evaluation of geometric phase, we need to know the dynamical phase. 
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From the formulas in ||^, one may find the relation 

Sniro) = -- / <n\H\n>dt = -- < n\{— + -Mw^x^)\n >dt 
n Jo n Jo 2,M I 

The geometric phase for the wave function in (22) under ro-evolution is therefore given as 

7 = vr[l + alD^ + ^ iBnl^n + (24) 

n=0 ^ O COS p 

If \Bn\ is given as |-B„| = 5„,m,, then one can find that the 7 given in (24) is equal to 
7m of (21) up to an additive constant 27i. 



5. Geometric phase under the half-period evolution 

Under a time evolution of ro/2, any solution of classical equation of motion of SHO changes 
sign with same magnitude. As noted in 0, the existence of quasiperiodic eigenfunctions 
depends on the periodicity of p(t) and Xp(t). Therefore, if we take Xp = for SHO (that is, 
in the absence of fictitious particular solution) , the geometric phases of the eigenfunctions 
can be defined under the evolution of half-period tq/2. Indeed, one may find that 

iJn{x, To/2 + t)= exp[-(l/2 + n)i7T]ij„{x, t). (25) 

The geometric phase for the eigenfunction ipn\xp=o under the half-period evolution is given 

as 

Ur»/2) = \-,Jr„) = („ + + 1±|1]. (26) 
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Due to the fact that phase is given up to an additive constant 27r, for the wave functions 

oo oo 
n=0 n=0 

geometric phases can also be defined under the evolution of the half-period. 

Making use of the previous formulas (particularly, (24)), one may easily evaluate the 
geometric phases of \E'e, ^o- In addition, \E'e(x, t), \E'o(x, t) has the following symmetry 

^e{-X, t) = ^eix, t) ^o{-X, t) = -^oix, t) . 

Due to the completeness, it may be stated that for a wave having the symmetry ^{—x, t) = 
±\E'(— x,t), the geometric phase can be defined under the half-period evolution. 



6. Perturbative Mathieu's oscillator 

In this section, we will consider the Mathieu's oscillator system described by the Hamil- 
tonian 

H = — + —{a + lQecos{2t))x\ (27) 

where M, a are positive numbers and the real number e will be assumed small. The 
classical equation of motion of the system is given as the Mathieu's equation 

x + {a + 16e cos(2t))x = 0. (28) 

The Mathieu's equation has long been studied in mathematics and, when a is close 
to n"^ {n = 0, 1,2, ■ ■ ■), some perturbative solutions are known. For example, when a = 
1 — 8e — 8e^ + 0{e^), a periodic solution (Mathieu function) of (28) is known as |TH 



oo 2,^(^^ 2^r+l^r+l 

cei{t, e) = cost + ^{7 — -r - 7 —7 — + 0{e'+^)} cos(2r + l)t. 

fr{\r + iy.r\ (r + l)!(r + l)! ^ ^ ^ 
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When a = 1 + 8e — 8e^ + 0{e^), there is another series solution sei(t, e) which reduced to 
sint in the hmit e ^ [15]. 

For the Mathieu's system, if u{t),v{t) are two hnearly independent solutions of (28), 
by adopting the definitions in (4), the eigenfunctions and their geometric phases are given 
as in (3) and (7), respectively For the calculation of geometric phase, instead of trying 
to find two linearly independent solutions, p(t) will be perturbatively evaluated. Similar 



method has been used in [T^ where, by examining the condition that the perturbative 
Lewis' invariant gives an ellipse in the phase space, the regions of (a, e) for the stable 
solutions of the Mathieu's equation are given in a reasonably good agreement with the 
exact results in the e — > limit. 
We rescale the p(t) as 

(29) 

Then one may find |l| that the p satisfies the following differential equation 

^+ (a + 16ecos(2t))p- ^ = 0. (30) 

The geometric phase for an eigenfunction ipn{x-,t) may then be written as 

ln={n + \) r^dt, (31) 
2 Jo 

if p is periodic under vr-evolution. 

We assume the perturbative expansion of p for small e as 

P = Po + epi + eV2 + e'p3 + 0(e^). (32) 
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From (30), one may find the differential equation for the po 

-.1 

Po + «Po - -33 = 0. (33) 
Po 

We will assume that po is constant, which gives 

Po = a-^/^ (34) 

The above assumption amounts to that we only consider the representation (s) of Mathieu's 
oscillator system which reduces to the stationary representation of SHO in the e — limit. 
Such assumption may be natural except for the cases a = (n)^ {n = 1,2, ■ ■ ■), since we 
should find, for the geometric phase, the representation of quasiperiodic eigenf unctions 
under the evolution of vr in the Mathieu's system (which reduces to the SHO in the e — > 
limit), while for r? the stationary representation is the unique one of quasiperiodic 
eigenfunctions under the vr-evolution in the SHO system of e = 0. 

By comparing terms of same order in (30), we find the recursive equations: 

cos 2t 

P;^+4api = — 16 — -7—, P2 + 40-/52 = —16pi cos 2t + 6a^''^p?, •••• (35) 
The solutions are given as 

4cos2t 

Pi = 1717 TT' 

a^i^\a — 1) 

4(5a-2) 4(5a-2) 
= aV4(a-l)^ + aV4(a-l)>-4) ^"^'^' ■■■■ ^'^^ 

In obtaining (36,37), we require that pi,p2 ■ ■ ■ are periodic functions of tt for general a. 

From (31), the geometric phase for the eigenfunction ipn in the representation mentioned 

above is given as 

_ 2 8(n+ |)7r 

v^(a-l)2+^^')- ^2^) 
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In evaluating (38), p2 has been used to show the absence of a term of the order of in 

In • 

The expressions of pi,p2 in (36,37) indicate that this perturbative method is not vahd 
when a = 1,4 (presumably, also when a = (n = 3,4- ■ ■)). In it has been proven 
that geometric phase can not be defined for a harmonic oscillator system if one of the 
classical solutions diverges as time goes to infinity. Indeed, it has been numerically proved 
that for the cases a = 1 or a = 4 with nonzero e, one of the solution diverges as time goes 
to infinite fl^ (see also |18|, H). The numerical study fl^ also indicates the existence 
of the representation (the existence of real characteristic exponent) where the geometric 
phase can be defined for small e and (see 0). 



7. Summary and discussions 

The general formalism for the geometric phases of harmonic oscillators ^ has been explic- 
itly applied for (driven) SHO. By calculating the geometric phases for the eigenfunctions, 
and relying on the completeness of a set of eigenfunctions, geometric phase for a general 
wave function of SHO is given. 

The explicit study on the SHO system suggests a new interesting feature on the ge- 
ometric phases of harmonic oscillator whose Hamiltonian is periodic with period tu\ 
For the case of two linearly independent periodic homogeneous solutions of period tu or 
2r/f, we can think of taking fictitious particular solution even for the undriven oscillator, 
which would give different geometric phases without changing the period needed for the 
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definition of geometric pliase. However, for the case where one of the homogeneous solu- 
tion is not periodic with period th or 2thi if we take fictitious particular solution for the 
F = oscillator system, the period needed for the definition of geometric phase should 
be changed or in some cases the geometric phase can not be defined. 

The geometric phase of an eigenfunction of SHO is interpreted as the contributions 
from the two different motions of probability distribution (pulsation of width and oscil- 
lation of the center). It has been shown in Q], the eigenfunction whose center oscillate 
is related to the coherent state constructed by applying displacement operator to the 



ground state of the stationary representation. As suggested in |[T^ , the eigenfunctions of 
harmonic oscillator is closely related to the generalized coherent states, and the eigenfunc- 
tions of probability distribution function whose width pulsates would be related to the 
generalized coherent states constructed by applying the squeezing operators to the ground 
state of the stationary representation |T3. The geometric phases has been studied for a 



generalized coherent states in It would be interesting to compare our results with 
those for coherent states. The contribution to the geometric phase from the motion of the 
center of probability distribution {xp) for the driven case is 0{l/h) in a general harmonic 
oscillator. For the F = case, the contribution from the oscillation of the center which 
is from taking fictitious particular solution is also 0{1/K). 

We also studied the geometric phase of perturbative Mathieu's oscillator system of 
small e. The representation in which we calculate the geometric phase must be the only one 
representation of the system where the geometric phase can be defined under vr-evolution, 
since the stationary representation is the only one representation of the corresponding 
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SHO (e = 0) system with quasiperiodic eigenfunctions except for the cases a 



(n = l,2, 



As a final remark, we note that geometric phase has been studied in the finite dimen- 



sional Hilbert space of SHO constructed from bra and ket vectors in the stationary 



representation, which is based on the Pegg-Barnett formalism of finite dimensional 



state space (see, also would be interesting to study the geometric phases for the 

finite dimensional Hilbert space in a representation other than the stationary one. 
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